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Task  I 

Modulation  of  Light  for  Optical  Generation  of  Sound 

by 

Henry  E.  Bass 

This  task  is  devoted  to  optical  generation  of  sound  at  audible  frequencies. 

The  optical  source  of  energy  roust  be  modulated  to  produce  a  time  dependent  heating 
and  this  time  dependent  heating  is  then  converted  to  a  pressure  variation.  We 
have  been  considering  conversion  from  light  to  sound  in  a  gas  (classical  spectro 
phone)  and  at  a  surface  (as  observed  by  Alexander  G.  Bell).  Our  efforts  during 
the  past  year  have  been  devoted  to: 

1.  Developing  the  theory  of  sound  conversion  at  a  surface,  preparing  sur¬ 
faces,  and  preliminary  measurements  to  compare  with  theory.  Our  results  to  date 
are  summarized  in  Appendix  A  which  follows. 

2.  Developing  a  theory  which  describes  transfer  of  energy  to  a  pressure  pulse 
in  a  gas  (  similar  to  the  liquid  problem)  and  comparison  of  predictions  to  experi¬ 
ment.  The  theory  is  described  in  Appendix  B;  experimental  results  have  not  been 
completely  analyzed. 

3.  Continuing  study  of  sound  amplification  in  a  reacting  gas. 

We  have  also  made  progress  in  the  study  of  pulse  density  modulation  as  a 
means  of  converting  a  string  of  rapid  pulses  to  an  acoustic  signal  (see  Appendix  C 


attached) . 
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APPENDIX  I- A 

Generation  of  Sound  at  a  Surface  Subjected  to  Modulated  Light 

Introduction 

Generation  of  acoustic  waves  by  absorption  of  electromagnetic  radiation 
indicent  on  a  surface  was  first  reported  by  Bell.1  More  recently,  this  tech- 
nique  has  been  adapted  to  studies  of  surface  absorption  coefficients.  The 
possibility  of  exploiting  this  phenomena  as  a  means  of  generating  ultrasound 
for  communications  purposes  was  raised  by  Westervelt  and  Larson. 3  They  treated 

the  case  where  thermal  conduction  can  be  ignored  which  is  true  at  very  high 
frequencies  in  water.  Their  result  yields  small  conversion  efficiencies 

except  at  very  high  frequencies.  Temkin4  treats  the  case  of  a  flat  surface 
with  its  temperature  moldulated  and  concludes  that  the  efficiency  for  conver¬ 
sion  of  evergy  into  sound  is  proportional  to  to  2  with  a  value  near  10~5  at  10Hz. 
This  would  imply  that  a  100  watt  light  bulb  turned  on  and  off  at  10  Hz  would 
provide  a  sound  pressure  level  of  125  dB  in  air.  This  is  not  the  case.  The 
reason  for  the  descrepancy  lies  in  the  neglect  of  the  energy  required  to  main¬ 
tain  an  average  energy  flow  away  from  the  surface. 

In  the  following,  we  will  treat  sound  generation  at  the  surface  of  a  plate 

coated  with  a  very  thin  layer  of  optically  absorbing  material  (lamp  black  or 
AIO2) •  We  will  be  interested  in  the  audible  frequency  range  particularly  the 
low  end  of  this  range  (M.0  Hz)  . 


Theory 

From  Reference  4,  writing  the  temperature  of  the  abosrbing  surface  as 


T  =  T  +  e  cos  to  t 
w  o 


(1) 


and  assuming  that  there  is  no  heat  flow  out  of  the  rear  plate,  the  acoustic 


intensity  becomes 
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I  =  %  p  0)K  C  B2£2  (2) 

o  o  o  o  v  ' 

where  <  is  the  thermal  diffusivity,  c  is  the  speed  of  sound,  and  (3  is  the 
coefficient  of  thermal  expansion  which  takes  on  the  value  1/T  for  an  ideal  gas. 

Burdick  and  Arnold  have  solved  for  the  temperature  fluctuations  in  a 
thin  film  absorbing  modulated  radiation  of  the  form 


I.(t)  =  I  /2<l-ejMt) 

X/  O 


(3) 


They  write  the  energy  balance  equation  in  the  form 

otl  =  C  p£(dT*/dt)  +  BT*  (4) 

P 

where  a  is  the  fraction  of  the  radiation  sbsorbed,  (C  p£)  is  the  specific 
heat  of  the  absorbing  surface,  T*  is  the  difference  between  the  surface  temp¬ 
erature  and  the  temperature  of  the  surroundings  and  B  is  the  effective  thermal 
conductance.  The  effective  thermal  conductance  is  the  power  transported  from 
the  surface  (unit  area  -°K)"1.  This  can  be  related  to  the  thermal  conductivity 
of  the  gas,  by 

B  -  k  /S  (5) 

o  o 

where  6^  is  the  distance  in  the  gas  necessary  to  establish  the  temperature 
difference  T*,  and  k  is  the  thermal  conductivity  of  the  gas. 


Eq.  (4)  has  a  solution 


T* (t)  = 


a!o  aIoexp(jmt) 
2B 


2(B+ja)C  p£) 


(6) 


For  the  purposes  of  Eq.  (2),  we  are  interested  only  in  the  real  part  of  the 


oscillating  term  in  Eq.  (6),  T*  (t) ,  so 

a.  c. 


T*  (t)  = 
a.c.  2 


-a  I  (B  cos  cut  +  wC  p£  sin  ut) 

_ o _ P _ 

[By-mz(C  pi)7] 


(7) 


We  see  that  in  Eq.  (.7),  there  are  two  oscillating  terms,  each  will  be  dominant 
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for  different  ranges  of  frequencies.  Prior  to  examining  the  intensity,  we  will 

establish  the  limiting  conditions  on  Eq.  (7). 

We  will  assume  6  is  the  thickness  of  the  thermal  layer  near  the  surface 
o 

% 

which  is  given  by  2tt(2k  /oj  )  so 


O  /  0)  x%  1/0)  ,  ^ 


B  =  (rf—)  2  =  f-  (  ~  p  C  k  p 

2tt  2k  2  tt  2  o  po  o 

o 


(8) 


If  we  let  p  C  =  2.83  x  10  ^  cal  (cm^K)*"1  and  k  =  6.22  x  10  5  cal  (cm3sec  K)  1 
o  po  o 


which  are  values  for  N2  at  300K, 

B  -  1.5  x  10”5  cal  (cm3K) (rn/sec) 


(9) 


The  quality  (C^p£)  is  more  difficult  to  estimate  since  our  surface  is  less 

well  understood.  We  will  adopt  for  gold  from  our  experiment  (see  next  section) 

which  gives  p£  -3.6  x  10“3  cal  (cm^K)“*.  Eq.  (7)  becomes 


—Oil  3*" 

/^x_  O  r1.5  X  10~5(a)/sec)  2cos  mt  +  3.6  x  10_3sinot1  ..-i,  ?TrN 

i*  „  Vw-  — r  -- ■  L - -  7T -= T'l h -TOT  "  - ~ - Jcal  1  ( cmz K)  (10) 

a.c.  2o)  2.25  x  10  lu/sec  -  o>  *2.3  x  10  b 


We  see  that  except  near  uj  =  0, 


T*  (t) 

a.c. 


al  sin  a>t 
o _ 

2a>  C  p£ 

P 


(ID 


The  sin  u>t  term  indicates  that  the  variation  in  gas  temperature  is  90°  out  of 
phase  with  the  light  source. 

We  can  now  rewrite  Eq.  (2)  as 


T  1  2j  2  POICQC0  1 

1  8  “  Z0  T2(C  pi)  2  u 

o  p 


(12) 


recognizing  that  when  0)  ->  0,  Eq.  (11)  is  no  longer  valid.  Writing  k  =k  /p  C  , 

000  po 


i  k  c  it 

T  _  1  2t  2  o  0  1  1 

1  8  “  To  C  T  7  (C  p£)^  to 

po  o  p 


(13) 


or,  assuming  plane  waves. 
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=  h  al  c  /T  (C 


(14) 


It  is  interesting  to  note  that  the  acoustic  intensity  goes  as  the  square 

of  the  optical  intensity.  Since  we  have  assumed  plane  waves,  this  means  that 

% 

the  acoustic  pressure  is  proportional  to  1^  and  inversely  proportional  to  w  . 

r\ 

This  result  can  now  be  compared  to  the  results  of  Roscenwaig  and  Gershoz 
and  Parker.6  Roscenwaig  and  Gersho  treated  the  more  complete  problem  of  a 
finite  backing  plate  and  an  absorbing  surface  of  finite  thickness.  If  we 
let  the  surface  thickness  approach  zero,  assume  all  the  energy  is  absorbed,  and 
the  plate  has  a  small  thermal  conductivity,  their  result  can  be  written  as 

«P(t)  »  QeK-t-V4)  (15) 

where 


YP, 


(31 


Q  = 


I- 


(1-j ) gk 


a  T  2k  (g2-a  2)  2k  oca  2/a 
ggo  s  s  gbsg 


] 


(16) 


Z  is  the  length  of  the  test  cell  (they  assume  a  finite  length),  g  is  the 
absorption  coefficient  which  we  assume  to  be  large,  kg  and  k^  are  the  thermal 

conductivities  of  the  surface  and  gas  respectively,  a^,  and  a  are  the  thermal 

y 

diffusivities  of  the  backing  plate,  surface,  and  gas  respectively,  a  =  (o)/2a  )  2, 

8  8 

y 

and  a  =  (1+j )  (o)/2as)  2.  Substituting  and  cancelling  terms, 


Q  = 


yP  g2I 
o  o 


(1-3) 


4£  T 


1  (o 2.  .  /  \  K  a,  a  /a  * 

g  o  (gz+  jrn/a  )  g  b  s  g 


k  a,  a  2 /  a  % 


0)2 


(17) 


Since  we  have  assumed  g2  is  large. 


(1-j) 


yP  I 
o  o 

y  y 

4£  T  k  oca  2/a  2  a) 2 
g  o  gb  s  g 


(18) 


Although  the  expression  is  quite  approximate  due  to  the  assumptions  made,  we 


note  that  6P  a  m  which  is  consistent  with  the  earlier  finding  that  I  a  aj~. 


Also,  6P  a  I  /T  consistent  with  our  earlier  result, 
o  o 
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For  a  fixed  length  cell,  Parker  finds,  for  the  pressure  midway  in  the  cell, 

(19) 


P 

o  o 


po 


Again,  6  is  proportional  to  a)  and  I  but  there  is  a  difference  between 
P  o 

Eq.  (20)  and  Eq.  (13)  in  the  way  that  C  enters.  Parker  finds  that  <5p  a  C  2; 

_> 

Eq.  (13)  gives  Sp  a  C  2.  This  descrepancy  arises  in  the  application  of  the 
boundary  condition  at  x  and  £. 

For  an  absorbing  liquid  ignoring  thermal  conduction,  Westervelt  and  Larson 
have  shown  that  the  total  acoustic  power  generated  is 


W  =  chi).B2P  2/16p  C  2 
o  o  po 

where  3  is  the  coefficient  of  thermal  expansion  and  a  is  the  attenuation 
coefficient  for  the  light  source.  Note  that  in  this  case,  W  is  proportional 


to  0)  . 
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Experimental 

Referring  to  Eq.  (13),  if  one  wishes  to  generate  sound  in  air  (our 

initial  goal),  the  parameters  which  can  be  controlled  are  a,  the  fraction  of 

absorbed  incident  radiation,  I  ,  the  incident  intensity,  C  p£,  the  specific 

o  p 

heat  of  the  absorbing  material,  and  a),  the  frequency  of  modulation.  Each  of 
these  parameters  have  been  explored. 

For  experimental  purposes,  the  surface  which  we  can  most  readily  control 

is  A10  and  AuO  .  Several  surfaces  were  prepared  with  both  A1  and  Au  vacuum 
x  x 

depostied  on  4  x  4  x  \  inch  glass  :plates.  The  amount  of  metal  was  controlled 
by  placing  different  amounts  of  metallic  wire  in  the  bell  jar.  By  placing  a 
shield  around  the  melting  wire,  an  estimate  of  the  fraction  of  metal  ending 
up  on  the  plate  was  determined.  In  addition,  the  surface  resistance  was 
measured  giving  another  indication  of  surface  thickness.  The  amount  of  oxida¬ 
tion  was  controlled  by  varying  the  pressure  in  the  bell  jar  when  the  metal 
was  heated  and  the  rate  at  which  the  metal  melted.  Our  goal  was  to  achieve 
maximum  absorption  at  the  surface  (incident  intensity/ transmitted  intensity) 

with  minimum  metal  (thereby  minimizing  C  p£) . 

P 

The  metal  wire  to  be  melted  was  attached  to  a  tungsten  filament.  The 
tungsten  filament  was  heated  by  controlling  the  current.  Several  problems  were 
encountered  which  made  the  resultant  coating  only  partially  reproducible. 

First,  if  any  metal  remained  on  the  filament  from  one  attempt  to  the  next, 
more  metal  was  available  for  coating  and  the  surface  would  sometimes  be  thicker 
When  the  tungsten  filament  was  heated,  the  pressure  inside  the  bell  jar  proved 
very  difficult  to  control  thus  giving  non-repeatable  levels  of  oxidation.  Heat 
ing  the  tungsten  filament  in  the  presence  of  oxygen  also  caused  the  filament 
to  oxidize  leading  to  a  very  short  lifetime.  All  these  problems  are  associated 
with  vacuum  deposition  in  a  reactive  atmosphere. 
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A  sufficiently  large  number  of  surfaces  were  prepared  to  give  us  a  reason¬ 
able  range  of  a  and  (C^p£)  .  The  conditions  for  the  gold  black  furface  are 
given  in  Table  I.  Each  piece  of  gold  wire  has  a  specific  heat  of  approximately 
2.3  joule/°K.  For  an  average  filament  to  plate  distance,  5/8  of  the  gold 
ended  on  the  plate  so  the  total  specific  heat  of  the  gold  deposited  was  about 
1-5  Joule/ °K.  The  plates  have  an  area  of  0.01m2.  However,  if  in  the  future 
we  wish  to  optimize  the  surface,  more  work  must  be  done  in  controlling  the 
deposition  process.  The  surfaces  and  deposition  conditions  are  given  in  Table  I. 

The  intensity  of  the  incident  radiation  and  frequency  were  varied  in  two 
ways.  First,  a  Gen  Rad  1531  AB  Stroboscope  was  used  as  a  light  source.  This 
source  delivers  a  pulse;  the  duration  and  energy  can  be  selected  from  three 
combinations:  0.8  ys,  6  x  105  candela;  1.2  ys,  3.5  x  106  candela;  or  3.0  ys, 

1.1  x  107  candela.  The  capacitors  in  the  stroboscope  emit  a  loud  click  when 
discharging  so  this  unit  was  mounted  in  a  box  with  a  plexiglass  window,  ,  The 
glass  plate  with  the  light  absorbing  surface  was  placed  in  a  second  box  with 
another  plexiglass  window  to  let  in  light.  A  B&K  %"  microphone  was  placed  in 
the  box  with  the  abosrbing  plate. 

Typical  waveforms  for  the  measured  acoustic  signal  are  given  as  Figure  1. 

It  is  now  of  interest  to  compare  our  experimental  observations  to  predic¬ 
tions.  First,  referring  to  Figure  1,  we  note  that  the  period  of  the  wave  is 
not  very  sensitive  to  the  length  of  the  incident  light  pulse.  This  resulting 
pulse  scales  in  pressure  approximately  linearly  with  the  total  power  of  the 
pulse,  in  accordance  with  Eq.  (14).  Substituting  estimated  values  for  the 
specific  heat  of  the  surface  (used  earlier  as  an  example)  and  assuming  a  funda¬ 
mental  frequency  of  5  kHz,  gives  p  «  1.0  x  1(T4  Pa.  Our  experimental  result 
gives  a  value  of  p  =  1.0  x  10“l  Pa.  We  are  not  now  in  a  position  to  determine 
the  origin  of  this  order  of  magnitude  discrepancy.  Assuming  the  theory  is  cor- 
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rect,  the  most  uncertain  quantity  in  Eq  *  (14)  is  (C^pjO  .  We  have  used  the 
value  for  gold  which  is  quite  likely  much  different  than  gold  oxide.  Further, 
we  assumed  that  5/8  of  the  gold  wire  found  the  plate.  This  fraction  could  be 
as  low  as  1/8.  Future  experiments  to  narrow  this  descrepancy  will  require  a 


more  careful  characterization  of  the  surface. 
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Conclusions 

At  this  time,  we  are  not  yet  in  a  position  to  determine  how  accurately 
Eq.  (14)  describes  the  sound  generation  process.  This  relation  will  allow  us 
to  reach  tentative  conclusions  concerning  the  viability  of  this  process  for 
generating  acoustic  energy.  Referring  to  Eq.  (14),  we  have  no  control  of 
V  p^,  or  CpQ*  This  leaves  1^,  the  incident  power,  (C^p£)  the  specific 
heat  of  the  surface,  and  to,  the  frequency. 

In  order  to  reconcile  experiment  and  theory,  we  have  to  assume  that 

C  pH  ^  1.1  x  10”4  cal  (cm  K)-1  which  would  imply  a  film  thickness  ^2y.  Our 
P 

experience  suggests  that  films  of  gold  black  less  than  ^  0.5y  will  become 
transparent.  This  means  that  C^p£  can  reduced  at  most  by  a  fraction  of  4 
using  geld  black.  Such  a  reduction  would  give  an  acoustic  pressure  of 
*  1.6  Pa  (98  dB)  at  5  kHz. 

The  frequency  dependence  in  Eq.  (1  )  suggests  that  this  method  of  genera¬ 
ting  sound  is  most  useful  at  low  frequencies.  By  decreasing  the  frequency  to 
50  Hz,  the  SPL  will  increase  by  20  dB.  This  result  is  quite  interesting  since 
conventional  sound  sources  tend  to  be  least  effective  at  low  frequencies. 

The  input  power  dictates  the  resultant  SPL  as  one  would  expect.  In  order 
to  generate  80  dB,  2.8  x  10 3  watts/m2  is  required.  To  generate  sound  over  an 
area  of  0.1  m^  280  watts  of  optical  energy  would  be  required.  This  is  an 
order  of  magnitude  higher  than  convential  sound  sources.  Assuming  we  can 
generate  an  absorbing  layer  of  gold  black  0.5u  thick,  at  50  Hz,  then  280  watts 
will  produce  a  SPL  of  118  dB  which  is  comparable  to  convential  sources.  At 
lower  frequencies,  opto-acoustic  conversion  should  be  more  efficient  than 
convential  mechanical  speakers. 

The  above  suggests  that  with  our  present  surfaces  opto-acoustic  generation 
of  sound  in  air  will  be  useful  only  when 
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a.  sound  needs  to  be  introduced  into  a  region  where  it  is  not  feasible 
to  place  mechanical  drivers  or  run  power  cords  (vacuum  systems,  low  EMI  test 
regions,  fiber  optic  communication  links,  etc.); 

b.  a  low  frequency  (^50  Hz)  sound  source  is  required  but  tuned  air 
horns  are  not  feasible. 

Our  major  effort,  then,  must  be  to  reduce  the  specific  heat  of  our  light 
absorbing  surface.  Before  doing  so,  we  plan  to  explore  Eq.  (14)  more  thorough¬ 
ly  to  determine  the  range  of  validity. 


Table  I 

B^B  Mi  B^B  BBB  ®Sara^ffiSist!iS^  ^BSrbiSf^laflS® 

(Gold  Black) 


Plate 

// 

Starting 

Pressure 

Amount 
of  Wire 

Transmission  of 

Surface 

Appearance  of 

Surface 

Comments 

//I 

50jj 

1  piece 

1  cm 

0.039g 

Very  transparent  coat 

About  1/2  transmission 

Gold-black  edges  with 
darker  spot  in  middle 

The  time  needed  to  vapo¬ 
rize  varies. 

n 

70-80y 

2  pieces 

1  cm  each 

0 .078g 

Very  little  transmission 

No  reflection 

Golden  brown  spot  in  middle, 
darker  black  edges 

in 

80y 

3  pieces 

2/3  cm  each 
0.078g 

Glass  almost  as  it  was 
before  coating 

Very  little  coating 

Filament  is  2  1/2"  higher 
from  plate  than  normal. 

#4 

8  Op 

3  pieces 

2/3  cm  each 

0 ,078g 

Center  is  clearer  than 
edges;  edges  high  trans¬ 
mission 

Edges  are  blue-black 

No  reflection 

Filament  is  1"  -  1  1/2  " 
higher  from  plate  than 
normal 

New  fi 

.lament  1/2  as  long  as 

filament  #1 

Golden  black  center  with 
streaks;  very  reflective 
except  streaks 

Filament  1  1/2"  off  plate 

in 

70-80y 

2  pieces 

1  cm  each 

0 . 078g 

Low  transmission 

#6 

70-80^ 

3  pieces 

1  cm  each 

0 .078g 

Transmission  is  low  where 
coat  is  thick;  not  as 
reflective  as  # 5 

Black  streaks  with  surface 
colors  varying  from  gold  to 
brown  to  blue  to  black 

Filament  was  2  1/2"  off 
plate 

#7 

8Cm 

2  pieces 

1  cm  each 

0 . 078g 

No  transmission,  especial¬ 
ly  at  center 

1st  coat  was  shiny  (reflec¬ 
tive)  surface;  2nd  coat 
thicker  nonref lective  coat 
very  black  center;  small 
gold  streaks  are  reflective 

Filament  2  1/2"  off  plate 

At  15V  for  3  min.  1  piece 
melted;  at  30-40V  for  1  mir 
other  piece  melted. 

Diameter  of  wire  -  0.0508  cm 


#2,  #7  are  the  "loudest"  plates 
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Table  II 

Characteristics  of  Stroboscope 


Strobe  Speed  Setting 

High 

Medium 

Low 

Pulse  Length  (ysec) 

0.8 

1.2 

3.0 

Luminous  Intensity  (Cd)* 

0.6  x 

106 

3.5  x  106 

11  x  106 

Power  of  Pulse  (total,  watts/m2) 

1.5  x 

102 

8.8  x  102 

2.8  x  103 

Energy/pulse  (J/m2) 

1.2  x 

10 -4 

1.0  x  10“3 

8.4  x  10-3 

Acoustic  pressure  (Pa) 

— 

22/-15.X  10“  3 

150/-110  x  10“3 

*Beam  Width  -  10°  =  0.17  Steradians 
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Figure  1 

Typical  Waveforms 


18  Nov  82 

Stroboscope:  690  RPM 

Rockland  Filter: 

Low  pass  80KHz 
High  pass  300  Hz 
Gain  20  dB 

Oscilloscope : 

Vertical  0.1  v/div 

Horizontal  0.2  msec/div 


10  Nov  82 
Stroboscope: 
Rockland  Filter 
Low  pass 
High  pass 
Gain 

Oscilloscope : 
Vertical 
Horizontal 


4170  RPM 

80KHz 
300  Hz 
20  dB 


50  fflv/ div 
0.2  msec/div 
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APPENDIX  I-B 

The  following  pages  represent  appendices  (draft  form)  to  a  paper  being 
prepared  on  spectrophone  measurements  of  vibrational  energy  transfer.  Please 
note  that  the  equation  symbols  and  appendix  letters  correspond  to  the  planned 
publication  independent  of  this  report. 
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Appendix  A 

Solution  of  Heat  Flow  Equations 


The  kinetic  equations  for  the  spectrophone  response  as 
time  have  been  solved  by  Bauer.1  He  assumed,  however,  that 
tion  to  the  cell  walls  is  very  slow  compared  to  the  rate  of 
deactivation.  In  this  case, 


a  function  of 
thermal  conduc- 
collisional 


P(t)  -  P(o)  =  (P/nTC  )  E(r.ur.x  1n  )  (1  -  e  w  vs) 


_-t/x 


A1 


where  r  is  the  i —  right  hand  eigenvector  of  the  matrix  [vkv(x  1  +  uu/ 

l  *^  ■*“  *^ 

00  2 

RCvTz)]  as  defined  by  Bauer,  u  is  a  vector  of  energy  states,  v  is  a  matrix 

of  stochiometric  coefficients,  x"*1  is  a  matrix  with  the  inverse  of  the 

00 

concentrations  in  the  diagonal,  is  the  non-relaxing  specific  heat,  T 

is  the  temperature,  n  is  a  vector  with  the  total  excitation  resulting 

from  the  input  laser  pulse  for  each  energy  state,  t  is  time,  and  x1  is 

vs 

the  relaxation  time  (eigenvalue)  associated  with  the  i—  eigenvector. 

This  expression  can  be  conveniently  rewritten  in  the  form 


P(t)  -  P(o)  =  (P/nTCv)Z6i(l- 


-t/T1  V 

e  vs) 


A2 


where  6^  is  a  relaxation  strength. 

For  purposes  of  the  experiment  described  here,  thermal  conductivity 
can  not  be  ignored.  The  equation  for  heat  flow  in  the  cell  can  be  written 


80 

at 


,82e  +  1  30.  ,  .  . 

k(37~  7a7)+0(r.t) 


A3 
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where  0  is  the  difference  between  the  gas  and  tube  temperatures  (the  tube 
temperature  is  assumed  to  be  constant)  .  We  will  assume  that  the  heat  gene¬ 
ration  term  0(r,t)  is  a  product  of  space  dependent  and  time  dependent  parts, 

0(r,t)  =  £  f (r)g(t) 

where  k  is  the  "thermal  dif f usivity"  and  k  is  the  thermal  conductivity 

(k  =  kC^p) .  The  function  g(t)  is  the  heat  generated  by  relaxation  processes 

which,  from  A2  is 

t/T* 

g(t)  =  E6 ,  (1  -  e  vs).  A5 

i  1 

Actually  Eq .  A5  is  more  exact  than  A2  since  the  derivation  leading  to  A2 

assumes  that  the  pressure  change  is  transmitted  to  all  parts  of  the  cell 

instantaneously  (acoustic  transmit  time  is  ignored). 

3 

Weulevesse  and  Genier  use  the  translational  part  of  the  thermal  conducti¬ 
vities  and  heat  capacities  in  Eq  A3  then  include  diffusion  of  vibrational  energy 
to  the  walls  by  adding  another  term.  We  will  choose  to  consider  diffusion 
of  vibrational  energy  to  the  walls  separately  since,  for  our  conditions,  this 
effect  is  almost  always  negligible.  Even  though  we  consider  diffusion  of 
vibrational  energy  separately,  we  still  use  the  equilibrium  values  of  k  and  c^ 
even  though  these  include  vibrational  contributions. 

To  a  good  approximation,  k  =  5/3  r\  (C^/m)  (Euken  Approx.)  which  gives 
k  =  5n/2Tnp.  This  equation  remains  true  without  regard  for  the  number  of  degrees 
of  freedon  (beyond  translation)  considered  so  the  use  of  equilibrium  values  of 
k  and  C  will  yield  the  same  results  in  the  important  exponential  term  which 

follows  as  a  more  rigorous  treatment  based  on  k  and  C 

trans  v,trans 

Prior  to  proceeding  further  with  the  solution  of  A3,  a  form  for  f  (r) 
must  be  assumed.  In  the  following,  we  will  assume  a  uniform  axial  excitation 
(we  have  already  assumed  uniform  excitation  longitudinally)  then  we  will  assume 
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a  Gaussian  distribution  across  the  cell  radius, 


Uniform  Excitation 


For  uniform  excitation,  f(r)  =  1  and 


o  ,  j  (r£  ) 

0  (r ,  t)  =  Z- 


ka  2  n[Ji(a^  )]2'  5  Jl(a?i2  &(t  )e  n  dt 

n  n 


A6 


where  5  is  a  root  of  the  equation  J  (£  a)  =0  which  arises  from  the  condi- 
n  on 

tion  that  0  =  0  when  r  equals  the  cell  radius  a  and  are  Bessel’s 

functions.  Evaluating  the  integral  over  time  gives 


r\  /  .  \  2k  ,  o  n-i  /  -t/x,  _K£2tN  o  ^n^ 

0(r,t)  =  7—  E.6.(k^t  -  1)  x(e  i  -  e  ) — —  •  y  r  r  . 

k  n,i  i  n  i  a£  Ji (a£  ) 

n  n 


A7 


We  now  assume  that  the.  acoustic  transit  time  is  small  compared  to  relaxation 
processes  so  the  pressure  observed  is  proportional  to  the  average  temperature 
rise  across  the  cell  radius  where 


0(t)  =  2^/^© (r ,  t)rdr/7Ta2. 


A8 


Integrating  gives 


4  K  i 

e(t)  -  -4  f  z. 


az  k  n,i  £  2 
n 


(k£2t.  -  1)  1  (e  t/,Ti  -  e 
n  l 


-Kf2 n *) 


A9 


which  is  the  desired  result.  Writing  out  the  pressure  dependent  terms 
explicitly , 


fl/.N  _  A  RT 

(t)‘^  v 


x.  RTk52T.  t/x 

I.  (— (e 


C  p 


-RTk£ 2  t/C  p 
-  e  ) 


A10 


or 


.  6;  C  p  -ETk{?t/C  p  -(- - RTk{2/C  p) C 

n 


T  . 


n  v 


-1)  All 
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It  can  be  seen  from  Eq.  All  that  the  effect  of  heat  conduction  is  to  make 
the  effect  deexcitation  rate  larger  and  pressure  dependent  and  make  the  relaxa¬ 
tion  strength  pressure  and  time  dependent.  We  can  see  that  even  clearer  if 

n  4  ^i  C  p 

we  write  RTk£^  =  1/t^  and  (k^T±-  -^p)_1  =  A±n  so  Eq.  All  becomes 

n 


-t/x"  -t(l/T  -  1/t  ) 

e(t)=n^iAin^  tC-1>* 


A12 


In  the  proceeding,  we  have  ignored  energy  transfer  to  the  cell  windows. 
Since  that  mechanism  will  always  increase  the  energy  loss,  we  would  expect  that 

will  be  less  than  (C  p/RTk^ 2)  .  We  have  also  ignored  radiative  transfer  of 
lc  v  n 

energy  to  the  cell  walls.  This  will  be  considered  next. 

Ideally,  heat  transfer  to  the  walls  fo  the  tube  would  be  included  by 
introducing  spontaneous  emission  into  the  relaxation  equations.  To  do  so 
introduces  an  inordinate  degree  of  difficulty  since  the  resultant  matrices 
loose  the  symmetry  which  made  an  analytical  solution  possible.  We  choose  to 
treat  the  effects  of  spontaneous  emission  less  rigorously  but  in  a  manner 
which  gives  an  analytical  solution. 

Spontaneous  emission  will  have  two  experimentally  observable  effects .  At 
high  pressure,  since  the  effects  of  thermal  conduction  decreases  with  increasing 
pressure,  spontaneous  emission  is  a  mechanism  for  transferring  the  heat  supplied 
to  the  gas  by  the  laser  to  the  cell  walls.  This  effect  will  parallel  thermal 
conduction  and  will  be  important  for  high  gas  pressures  and  long  times.  At 
low  pressures,  spontaneous  emission  will  offer  a  parallel  relaxation  path  to 
the  slow  collisional  processes. 

On  a  time  scale  much  longer  than  required  for  collisional  energy  transfer, 
the  vibrational  modes  are  in  equilibrium  and  can  be  described  by  a  single  temp¬ 
erature.  In  the  absence  of  conduction  or  radiation  losses,  this  temperature  is 
given  as 
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T  =  T  +  (P/nTC  °°)Er.ur .x-1n  A13 

gas  v  l ~  -ex 


where  T  is  the  temperature  of  the  cell  walls. 

The  rate  of  heat  transfer  per  unit  length  to  the  walls  of  a  cylinder  of 

4 

gas  can  be  written  approximately  as 


Q  =  47TOa  T30R 
g 


A14 


where  a  is  the  Stefan-Boltzman  constant,  is  the  absorption  coefficient  for 
the  gas,  0  is  the  difference  between  the  gas  and  wall  temperatures,  and  R  is 
the  cell  radius.  We  have  assumed  that  the  absorption  coefficient  equals  the 
emissivity  and  have  ignored  terms  in  or  greater. 

The  rate  of  change  in  temperature  is  given  by 
-47TOa  T3R*£ 


e  =  -^r 


e 


A15 


where  i  is  the  cell  length  and  is  the  total  heat  capacity.  In  terms  of  the 
specific  heat/unit  volume, 

aa  T30 


0  =  - 


C  R 
v 


A16 


which  has  a  solution 

-th 

0  =  Ae 


A17 


where  x  =  C  R/aa  T3.  Note  that  x  is  no  more  than  the  lifetime  of  temperature 
r  v  g  r 

displacement  from  equilibrium  which,  at  high  pressure,  will  be  the  lifetime  of 

of  the  shortest  lived  vibrational  state.  Also  note  that  both  C  and  a  increase 

v  g 

linearly  with  the  pressure  in  the  pressure  broadened  regime  so  x^  will  be  inde¬ 
pendent  of  pressure. 

At  very  low  pressures,  Eq.  (A17)  will  fail  for  several  reasons.  First,  we 
can  not  then  assume  that  all  vibrational  modes  are  in  equilibrium  with  transla- 
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tion  .  Second,  the  gas  will  become  optically  thin  which  will  make  pressure 
dependent.  Fortunately,  at  low  pressures,  thermal  conduction  will  generally 
dominate  the  long  term  decay  and  Eq.  A17  will  represent  only  a  second  order 
effect . 

When  the  pressure  is  very  low  collisional  processes  are  comparable  in 
speed  to  spontaneous  emission.  As  a  result,  the  observed  relaxation  times  will 
be  a  combination  of  collisional,  spontaneous  emission,  and  wall  deexcitation 
(discussed  later) .  The  combined  rate  of  deexcitation  will  then  be 

3  =  K  +  K  +  7T 
c  w  p 

Where  Kj,  is  the  collisional  deactivation  rate,  is  the  rate  of  deexcitation  by 
the  walls,  and  is  the  rate  of  deactivation  by  spontaneous  emission.  There 
will  be  a  different  3  for  each  energy  level;  for  a  multi-level  system,  Eq.(18) 
is  only  useful  if  tt^  were  the  same  for  each  level  (it  is  not). 

To  a  first  approximation,  we  can  identify  the  relaxation  times  which 
result  from  the  general  reaction  matrix  (see  Eq.  A2)  with  specific  levels.  For 
example,  the  relaxation  process  which  contribures  most  to  sound  absorption  is 
associated,  to  a  large  extent,  with  the  lower  vibrational  energy  state  (assuming 
series  relaxation) .  The  relaxation  process  which  contributes  most  to  fluores¬ 
cence  decay  is  associated,  to  a  large  extent,  with  the  IR  mode  which  fluoresces. 
These  assignments  are  not  rigorous  and,  should  only  be  used  when  consider¬ 
ing  second  order  effects. 

With  this  assignment,  we  will  write 

1A  ,eff  =  1/T.  +  1/T.rad 
i  ii 

where  is  the  collisional  relaxation  time  in  the  absence  of  spontaneous  emission, 
rad  ,  , 

T-j_  is  the  radiative  lifetime  of  the  level  associated  with  t^.  Then  we  replace 
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i  ef  f 

in  Eq .  (A2)  with  .  The  energy  which  goes  into  spontaneous  emission  does 

not  cause  the  gas  temperatures  to  change  so  6^  is  replaced  by 


x  p  • 

There  are  a  variety  of  other  heat  loss  mechanisms  that  will  help  the 
gas  return  to  equilibrium  on  the  long  time  scale.  Some  of  these  include  heat 
conduction  to  the  windows  and  absorption  of  propagating  acoustic  modes.  The 
significant  factor  is  that  these  loss  mechanisms  will  increase  as  l/p(heat  con¬ 
duction)  or  be  independent  of  pressure  (acoustic  loss  upon  reflection)  so  they 
can  be  included  in  x  or  as  appropriate.  Our  method,  then,  will  be  to 
modify  x  as  necessary  to  give  the  correct  long  term  decay  at  low  pressure, 
and  t  to  give  the  proper  long  term  decay  at  high  pressures.  These  second 
order  corrections  should  provide  only  third  order  effects  on  the  more  rapid 
part  of  the  waveform  at  intermediate  pressures. 
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Appendix  B 

Effect  of  Wall  Collisions 

£ 

According  to  Margottin-Machou ,  Dogonette,  and  Henry,  the  measured 

relaxation  rate  g(=  1/t  )  can  be  written  as 

exp 

6  =  K  +  K  +  7T  B1 

c  w  p 

where  K  is  the  collisional  deactivation  rate,  K  is  the  rate  of  deexcitation 

w 

by  the  walls,  and  tt^  is  the  rate  of  deactivation  by  spontaneous  emission 
(including  trapping)  .  They  further  show  that  for  a  long  narrow  tube, 

K  =  y2(D/R2)  B2 

W 

where  D  is  the  diffusion  coefficient,  R  is  the  cell  radius  and  y  is  a  solution 
of 


pj1 (y)  =  hJ  (y) 
o 

where 

h  -  (Rv/2D) (a/2  -  a) 

where  v  is  the  mean  speed  of  molecules,  and  a  is  the  probability  of  de¬ 
excitation  during  a  collision  with  the  cell  walls  (^0.4  for  C02) .  The  ratio 
R/h  is  the  minimum  distance  a  molecule  must  travel  close  to  a  wall  before  being 
deexcited  —  A/a  where  A  is  the  main  free  path*  The  diffusion  coefficient  D  is 
pressure  dependent  so  we  will  write  =  DP. 

Margottin-Machou,  Dogonette,  and  Henry  found  that  for  pressures  of  1 
torr  and  above,  y  approaches  a  limiting  value  of  2.4  so  Eq.  B2  becomes 
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K  =  (2.4)  ZD  /.R2p 
w  o 


B3 


where  D  is  the  kinetic  self-diffusion  coefficient* 
o 

The  interesting  feature  of  the  derivation  leading  to  Eq.  B3  is  the 
fact  that  all  excited  states  are  deexcited  during  wall  collisions  with 
about  the  same  probability  so  we  can  write  (ignore  spontaneous  emission) 


1/ t 


i 

exp 


1/t 


i 

coll 


+  K 
w 


where  the  i  denotes  different  relaxation  processes.  We  can,  therefore,  sub¬ 
stitute  into  Eq .  A12  1/t^^  wherever  1/t^  appears  and  correct  for  the  change 
in  resulting  from  wall  collisions.  For  the  spectrophone ,  however,  a  change 
in  t  is  not  the  only  effect  on  the  response.  That  energy  transferred  to  the 
tube  walls  will  not  contribute  to  the  temperature  and  pressure  rise  in  the 
cell  so  the  magnitude  of  the  response  will  change  along  with  the  time  constant. 

The  energy  which  is  lost  to  wall  collisions  is  K  /  (K  +  K  )  times  the  total 

c  c  w 

energy  transferred  from  vibration  so  K  /  (K  +  K  )  should  give  the  fraction  of 

c  c  w 

the  total  energy  which  is  transferred  to  translation  due  to  collisional  pro¬ 
cesses.  We  must,  therefore,  replace  <5.  in  Eq.  All  with  6  (K  /K  +  K  ). 

i  x  c  c  w 
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Appendix  C  -  Acoustic  Effects 


There  are  several  purely  acoustic  effects  which  will  affect  the  spec- 
trophone  response.  These  include  "acoustic  transit",  excitation  of  normal 
modes  of  the  cavity,  and  propagating  waves.  We  will  consider  a  cylindri¬ 
cal  tube  as  shown  in  the  following  diagram. 

Z 

"m  1 micro phone 

1  * 


First  consider  effects  of  acoustic  transit.  If  the  tube  is  uniformly 
irradiated  in  the  r  direction,  the  temperature  at  all  radial  distances  for 
a  given  z  will  evolve  together.  The  transducer  only  senses  pressure  in¬ 
creases  at  r  =  a  but,  since  the  gas  is  not  allowed  to  expand,  this  will  be 
the  same  as  the  pressure  elsewhere  along  the  cross-sections.  There  will 
be  no  radial  wave  propagation  and  the  microphone  will  accurately  repro¬ 
duce  the  pressure  response.  Uniform  excitation  is  experimentally  diffi¬ 
cult  to  achieve.  As  a  worse  case,  consider  a  Gaussian  excitation  with  no 
irradiation  near  the  tube  walls.  In  this  case,  the  center  will  recieve 
the  greatest  energy  and  a  radial  pressure  wave  will  move  outward  impinging 
upon  the  transducer  a  time  t  =  a/C  following  excitation.  One  could  inte¬ 
grate  across  the  tube  cross-section  to  find  the  resultant  pressure  response 
but  such  an  exercise  would  only  be  useful  if  I(r)  was  known.  It  is  obvious, 
however,  that  for  the  Gaussian  case,  the  measured  pressure  will  contain  in¬ 
formation  about  the  kinetics  which  is  an  average  over  events  happening  on 
a  time  scale  of  a/C.  We  will  call  this  time  the  acoustic  transit  time  x 

a.  t 

Note  that  this  represents  a  worst  case.  In  practice,  the  beam  is  probably 

more  nearly  uniform  than  Gaussian  so  x  <  a/C. 

a .  t 
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There  is  a  second  effect  of  non-uniform  excitation  across  the  tube, 
that  is  the  excitation  of  radial  modes.  The  cutoff  frequency  for  the 
first  unsymmetric  mode  is  approximately  f  =  0.586  C/a.  Radial  modes  will 
enhance  some  components  of  the  pressure  response  and  diminish  others.  By 
setting  2?rf^  =  1/x  ,  we  can  represent  the  upper  frequency  response  by  a 
time  constant  x^  which  is  the  shortest  time  period  for  which  radial  modes 
will  not  significantly  impact  our  results.  Solving  for  t  , 

xc  =  a/1.172TTC.  Cl 

which  is  always  less  than  x 

a .  t 

The  above  suggests  that  for  times  less  than  a/C,  the  uniformity  of 
the  excitation  beam  can  begin  to  have  an  effect  on  the  measured  pressure 
response.  For  time  segments  greater  than  a/C,  acoustic  transit  and  radial 
modes  will  not  affect  the  resolution  of  the  experiment.  We  will  choose  a 
such  that  x  is  always  smaller  than  the  resolution  desired.  The  alter- 

3  •  t 

native  would  be  to  insure  a  uniform  beam  cross-section  out  to  the  tube  walls. 

As  the  laser  beam  propagates  down  the  tube  it  will  be  absorbed  in 
accordance  with  Beer's  Law, 

T  T  ~aZ 

1  =  V  ’  C2 

where  a  is  an  extinction  length.  There  is  no  way,  experimentally,  to 
remove  this  Z  dependence  of  excitation.  We  can  force  a  +  0  by  adjusting 
the  laser  frequency  and/or  gas  pressure  but  this  reduces  the  observed 
signal  proportionately  (to  first  order). 

As  a  result  of  the  exponential  illumination,  more  energy  will  be 
absorbed  near  the  entrance  window  giving  a  greater  pressure  rise  there 
than  elsewhere  down  the  tube.  This  difference  in  pressure  rise  down  the 
tube  will  cause  a  pressure  wave  propagating  down  the  tube  so  that  the 
signal  received  at  the  microphone  will  be  a  sum  of  signals.  Further,  if 


30 


the  tube  is  terminated  at  the  opposite  end,  this  pressure  wave  will 
bounce  off  the  opposite  end  and  set  up  standing  waves.  These  standing 
waves  can,  however,  be, effectively  eliminated  by  making  the  tube  long 
enough  so  that  L/C  is  much  longer  than  the  time  of  interest  for  studying 
kinetics.  Alternatively,  an  absorber  could  be  placed  at  the  exit  end  of 
the  tube. 

Assuming  all  pressure  changes  are  small  so  that  disturbance  propagated 
with  the  speed  of  sound,  the  pressure  received  at  any  point  z  will  be 


p(t,z) 


z  9p  (t ' , z 1 ) dz’ 
o  9z' 


C3 


where  t'  =  t  -  (z  -  z')/C  is  the  retarded  time.  Since  p(t',z')  is  directly 
proportional  to  the  energy  absorbed  (proportionality  constant  is  time  de¬ 
pendent  but  not  z  dependent),  3p(t ' ,z' ) /9z '  =  f(t')al  e_az  so 


p(t,z)  =  al  /Zf(t')e  az,dz'. 
o  o 


C4 


Further,  since  the  microphone  has  a  finite  length,  £,  along  z,  the  pressure 
must  be  averaged  along  the  microphone  length  so 

P(t)  =  j/Zm  +  £  p(t,z)dz  i 

Equations  C4  and  C5  can  be  integrated  to  give  p(t)  or,  z  and  £  can  be 

m 

selected  to  minimize  z  dependence.  This  latter  requires  that  z  -*■  0  and 

m 

1  %  T_  t  c •  The  first  constraint  is  readily  met,  the  second  is  not  since 


a  very  small  microphone  will  reduce  sensitivity. 
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Appendix  D 

Response  of  Microphone  and  Electronics 

In  addition  to  thermal  conduction,  energy  transfer  to  the  walls  of 
the  cell,  and  acoustic  resonances  of  the  cavity,  the  response  of  the 
microphone  and  electronics  can  affect  the  recorded  pressure  response. 

The  microphone  response,  at  low  frequencies,  is  determined  by  leakage 
of  air  from  front  to  back  of  the  diaphragm.  At  high  frequencies,  the 
diaphragm  begins  to  oscillate  in  a  resonant  mode.  The  frequency  re¬ 
sponse  of  the  electronics  can  be  controlled  by  equipment  selection  and 
variable  frequency  filters.  It  may  be  more  convenient  to  control  the 
upper  and  lower  frequency  cutoffs  with  filters  of  known  frequency  re¬ 
sponse  rather  than  to  rely  on  theoretical  predictions  of  the  frequency 
response  of  individual  elements.  This  is  especially  truewhen  the  micro¬ 
phone  or  other  components  might  introduce  frequency  dependent  phase 
shifts . 

The  voltage  output  of  a  simple  capacitor  microphone  can  be  written 

as 

v( t)  =  M(pc  -  Pb) 

where  v(t)  is  the  time  varying  output  of  the  microphone,  p  is  the  pres- 
sure  in  the  cavity,  p^  is  the  pressure  behind  the  diaphragm  and  M  is  a 
sensitivity  coefficient  which  should  be  independent  of  frequency  at  fre¬ 
quencies  well  below  the  microphone  resonance.  The  pressure  behind  the 
diaphragm  will  equilibrate  with  the  front  pressure  with  some  time  con¬ 
stant  dictated  by  the  resistance  to  air  flow  front  to  back.  We  can 
write,  then 

dp  /dt  =  (p  -  P,)/t 
b  cbm 


D1 


D2 
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and  rewrite  Eq .  D1  as 


(p  -  v/M)  =  v/Mt 
at  c  m 


D3 


or 


—  +  —  =  M  dp  /d t . 

dt  x  *c 

m 

The  function  p ^(t) ,  the  pressure  response,  can  be  determined  from  the 
preceeding  theory.  M  is  just  a  constant  which  can  be  determined  (if 
desired)  from  microphone  calibration  curves.  We  will  write  the  quantity 
M  dp^/dt  =  f(t),  then 

v(0  =  e~t/Tlu/’t  f(t,)et'/Tmat'. 

o 


D4 


D5 


Next,  we  can  envision  applying  this  signal  to  the  input  of  our 
electronics.  The  following  circuit  diagram  represents  this  situation. 

+  i 


v.  o- 

in 


+  1 


-o  v 


out 


R 


The  values  of  C  and  R  are  taken  from  manufacturers  specifications. 

equation  for  voltage  becomes 

v .  =  q/C  +  Ri 

m 

where  q  is  the  charge  on  the  input  capacitor  and  i  is  the  current. 

Recognizing  that  i  =  dq/dt  =  v  /R,  D6  becomes 

out 


The 


D6 


dv.  /dt  =  1/C  (v  /R)  +  dv  /dt 
in  out  out 


D6 


or 


dv.  /dt  =  dv.  /dt  +  v  /t  , 
in  out  out  e 


D7 


where  x  =  RC . 
e 
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Combining  Eqs.  D4  and  D7, 
dv 


out 

dt  *  vout' 1 e  '  "in' "m 


+  v  /t  +  v.  /t  =  M  dp  /dt 

mif-  Ci  -fm  m  r  /- 


D8 


Note  that  if  *i  <<  a  condition  which  can  be  readily  controlled  experimentally. 


dv 


+  1/t  v  =  M  dp  /dt  =  f(t) 
dt  e  out  c 


D9 


which  has  a  solution 


vout  (t>  ■  e“t/Te 


DIO 


Ignoring  the  effect  of  acoustic  modes,  f(t)  =  0(t)  from  Eq.  A12  so 

:tc)  -  l)e+t’y  ledt’ .  Dll 


v  (t)  =  e 
out 


"tAe  ft  ,  A  -l/xjj  _  +t'/xe^, 


/  Z.A.  e  '  tc (e 
o  n,i  xn 


Integrating,  Dll  becomes 


-t/T, 


T.  -  T 
,  1 


=  e  ,_[(-* - — )(1  -  e  + 


out 


_  -t(l/Ti  -  1/Te). 


n,  i  m  t  -  t  .  ‘ 
e  i 


.  TtcTe)(l  -  tc  ~  l/Te)')-j^ 

T  Tn 

e  tc 

t  /  T 

Now,  if  «  Tfcc>  the  second  term  becomes  -  (1  -  e  e)  and  the  effects 

of  thermal  conduction  are  effectively  filtered  out.  Although  we  can  typi- 

ally  meet  the  approximation  necessary  for  D9,  at  low  pressures  (•v  10  torr) 

Tt  so  to  make  T  «  T  would  make  t  «  t.  which  would  eliminate  t 

-L  e  tc  cl  i 

from  Eq .  D12,  all  kinetic  information  would  be  lost.  At  higher  pressures, 
Te  <<:  Ttc  a  reasonable  approximation. 


i 
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APPENDIX  I-C 

Pulse  Density  Modulation  of  Light  for  Optical  Generation  of  Sound 

by 

Henry  E.  Bass 

Sound  can  be  generated  at  a  surface,  in  a  gas,  or  in  a  fluid  when  optical 
energy  is  converted  to  a  pressure  rise.  The  mechanism  for  this  conversion 
can  take  several  forms  ( thermoelastic  expansion,  dissociation,  vaporization, 
etc.)  but  we  will  be  interested,  primarily,  in  thermal  expansion  which  is  domi¬ 
nant  in  a  liquid.  The  thermal  expansion  mechanism  gives  a  pressure  rise  which 
is  proportional  to  ~  (3V/3T)  .  The  ~  term  gives  rise  to  a  1/T  term  where  T  is 

V  S  V 

the  pulse  duration.  The  1/T  dependence  leads  us  to  conclude  that  maximum  sound 
will  be  generated  by  a  series  of  short  pulses.  If  we  are  to  make  the  resultant 
meaningful,  these  pulses  must  be  spaced  in  occurance,  with  minimum  length,  and 
an  amplitude  such  that  an  audible  signal  is  produced.  We  could  amplitude  modu¬ 
late  a  string  of  pulses  to  achieve  this  goal  but,  from  a  practical  point,  it  is 
much  easier  to  modulate  the  rate  of  pulses  than  it  is  their  amplitude.  For  that 
reason,  we  have  chosen  to  examine,  in  detail,  pulse  density  modulation.  Pulse 
density  modulation  spaces  a  sequence  of  pulses  with  a  time  between  pulses  pro¬ 
portional  to  the  input  signal  amplitude  (to  some  power) . 

If  we  assume  a  rectangular  pulse,  we  can  control  the  time  between  pulses 

and  pulse  length.  We  want  to  produce  an  acoustic  output  which  will  be  a  true 
representation  of  the  input  electrical  signal  but  it  is  not  obvious  how  one 
can  best  accomplish  this  objective.  For  discussion  purposes,  let  us  assume 
that  we  wish  to  reproduce  a  sine  wave.  Further,  as  a  starting  point,  we  wall 
assume  that  the  optimum  output  is  obtained  when  the  number  of  pulses/unit  time 
is  proportional  to  the  input  signal  amplitude  at  that  same  point  in  time,  i.e.. 


Np  (t)  -  C  V(t), 


(1) 


36 


where  Np(t)  is  the  number  of  pulses/unit  time,  and  V(t)  is  the  input  waveform. 
Eq .  (1)  is  predicted  on  the  assumption  that  the  output  waveform  will  vary 
linearly  with  energy  input  to  the  surface  provided  the  input  to  the  surface 
consists  of  a  number  of  pulses  of  the  same  length. 

For  a  sine  wave  output,  V(t)  =  sin  cot,  so 

Np ( t)  =  C  sin  cot.  (2) 

The  number  of  the  pulse  can  never  be  negative  so  we  must  modulate  about  some 
average  number  of  pulses.  A, 

Np ( t)  =  “  (a  +  b  sin  cot).  (3) 

Note  that  although  a  +  b  =  2,  the  ratio  a/b  controls  the  degree  of  modulation. 
With  a/b  =  1,  there  are  no  pulses  emitted  when  V(t)  has  its  most  negative 
value;  when  a/b  =  0,  there  is  no  modualtion.  We  will  require  that  a  ^  b  so  the 
term  in  parenthesis  will  not  go  negative. 

Expressed  in  terms  of  the  time  between  pulses,  t,  , 

b 

2 

=  —  (  a  +  b  sin  cot)"*1.  (4) 

Let  Pr  =  maximum  pulse  repetition  frequency,  then 

4 

%  =  (a  +  b  sin  wt)"1  (5) 

if  we  want  a  maximum  duty  cycle  of  50%.  If  we  want  a  smaller  duty  cycle,  f^, 
^(t)  =  a/ff^^Ca  +  b  sin  wt)]  (6) 

where  f  is  the  maximum  fraction  of  time  the  light  source  is  on. 

In  order  to  examine  the  waveform  generated  by  a  series  of  pulses  spaced 
in  time  according  to  Eq  (6),  a  numerical  simulation  was  employed.  We  examined 
one  cycle  of  a  100  Hz  sine  wave  (V(t)  =  sin  2  Tift).  We  further  assumed  that 
the  pulse  duration  was  10  ysec  or  10-3  cycles  (actually,  we  used  1024  to 
simplify  the  Fourier  Transform)  .  We  then  place  a  1  or  a  0  in  a  1024  point 
array  to  indicate  a  pulse  on  or  off.  We  did  this  by  evaluating  Equation  (6) 
and  determining  if  the  time  since  the  last  1  was  greater  than  or  equal  to 
tb(t)  for  each  10  y  sec  interval  of  time.  This  array  was  then  transformed  into 
the  frequency  domain  to  determine  the  resultant  waveform. 
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There  were  two  features  of  the  frequency  domain  spectrum  of  interest; 
the  signal  to  noise  ratio  and  the  fraction  of  energy  appearing  at  the  desired 
frequency  (100  Hz) .  We  took  the  signal  to  noise  ratio  to  be 

20  log! o [A(100Hz) /A(200Hz) ]  where  A(f)  is  the  amplitude  of  the  Fourier  Transform 
at  a  frequency  f.  We  took  to  conversion  efficiency  (CE)  to  be 
20  log ! o [ A(100Hz) /A (DC) ]  since  most  energy  lost  to  the  100Hz  signal  appeared 
in  the  DC  term. 

The  first  set  of  results  (Figure  1)  are  for  a  pulse  length  of  1/6400  sec. 

It  can  be  seen  that  the  S/N  ratio  is  reasonably  constant  for  b/a  ^  0.5  to  1.0 
but  it  has  a  low  value  ^10  dB.  The  conversion  is  also  given  for  a  pulse  length 
of  1/12,800  sec;  the  S/N  is  better,  the  CE  is  worse.  Note  that  for  these  figures, 
the  duty  cycle  has  a  maximum  of  50%. 

The  next  results  (Figure  2)  are  for  a  pulse  length  of  10  psec  but  variable 
duty  cycle  fraction.  First,  looking  at  the  S/N,  we  see  that,  with  the  exception 
of  some  peaks,  the  S/N  varies  directly  with  PRF  from  about  22dB  for  a  PRF  of 
6400  Hz  (64  x  the  fundamental  frequency)  to  'v  10dB  for  a  PRF  of  1600  Hz  (16  x 
the  fundamental  frequency) .  The  conversion  efficiency  does  not  change  much  in 
this  PRF  range;  it  drops  uniformly  with  the  b/a  ratio  (at  b/a  =0,  CE  =  0) . 

The  plots  of  S/N  do  illustrate  that  optimum  S/N  occurs  for  b/a  <  1.0.  This  means 
that  when  the  signal  swings  negative,  the  number  of  pulses/sec  emitted  by  the 
light  source  never  goes  to  zero.  The  optimum  swing  in  signal  (b/a)  is  dependent 
on  PRF;  it  decreases  as  the  PRF  increases. 

The  above  results  suggest  that  PDM  will  be  effective  (acceptable  S/N)  only 
if  the  pulse  length  is  short  and  the  PRF  is  high.  Considering  light  sources  off 
the  shelf,  this  argues  for  a  fast  Xenon  flash  lamp  or  a  laser  diode  system.  A 
typical  off  the  shelf  Xenon  flashtube  emits  a  maximum  energy/flash  of  15  joules  in 
<  15  ysec  at  a  maximum  of  2500  pps.  A  laser  diode  system  will  emit  an  average 
of  6  m  watts  with  a  pulse  length  less  than  10  nanoseconds  at  100  M  pps. 


Puise  Le/ig+h  -  s/(t>Hoc>  sc  c 
Max  P  RF  ~  I &00  Hz. 
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Task  II 

This  task  involves  the  generation  of  acoustic  waves  from  the  interaction 
of  a  modulated  CO2  laser  pulse  with  a  highly  absorbing  liquid  medium  such  as  water. 
Earlier  experiments  in  water  with  proton  beams  [Hunter  et.  al.,  J.  Acoust.  Soc. 

Am.  69 ,  1557-1562  (1981)]  and  with  laser  pulses  [Hunter,  Jones  and  Malbrough,  J. 
Acoust.  Soc.  Am.  69,  1563-1567  (1981)]  gave  both  dipolar  and  tripolar  signals 
that  did  not  vanish  at  4°C. 

Since  these  pulses  were  presumed  to  be  generated  by  the  therraoelastic  properties 
of  the  liquid,  it  was  not  understood  why  they  still  existed  at  4°C.  Recently,  Pierce 
[J.  Acoust.  Soc.  Am.  _72,  S13  (1982)]  has  suggested  that  acoustic  signals  may  be 
generated  at  4°C  if  one  considers  the  effect  of  volume  relaxation.  He  has  suggested 
a  series  of  experiments  to  test  these  predictions. 

We  have  increased  the  peak  output  of  our  CO2  laser  system  with  the  acquisition 
of  an  optical  scanner  that  will  enable  us  to  vary  the  prf  from  dc  to  2  kHz.  We 
hope  soon  to  obtain  data  on  the  amplitude  of  the  thermoelastic  signal  as  a  function 
of  temperature  and  a  variety  of  other  liquid  parameters. 


Task  III 
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Propagation  of  Sound  through  Gases  with  ;m  Overpopulation  of  Excited 

Vibrational  States. 

Introduction 

The  propagation  of  sound  through  gases  under  a  variety  of  non-equi¬ 
librium  conditions  has  been  studied  the  past  ten  years.  These  studies 
have  predicted  sound  amplification  in  the  presence  of  a  chemical  reac¬ 
tion1  4  or  optical  pumping  of  vibrational  states5  7.  Bauer  and  Bass 
have  named  the  later  process  "SACER",  which  is  an  acronym  for  "sound 
amplification  from  controlled  excitation  reactions."  The  predicted 
amplification  is  attributed  to  the  preferential  heating  of  the  compres¬ 
sion  part  of  the  sound  wave  either  by  the  pumping  process  or  the  relaxa¬ 
tion  process  or  both.  The  amplification  due  to  pumping  results  from  the 
selective  absorption  of  pumping  energy  in  the  compression  part  of  the 
sound  cycle.  The  amplification  due  to  relaxation  results  from  the  tem¬ 
perature  dependence  of  the  vibrational  energy  transition  rate.  The 
increase  in  the  rate  of  v-t  energy  transfer  with  temperature  causes  the 
relaxing  vibrational  energy  to  be  selectively  "dumped"  into  the  hot  (or 
compressional)  part  of  the  sound  cycle.  When  the  sound  wave  is  propa¬ 
gated  in  the  presence  of  pumping,  both  effects  are  likely  to  be  opera¬ 
tive  and  the  two  will  be  difficult  to  separate.  In  addition,  the  non¬ 
equilibrium  conditions  produced  by  the  pumping  make  accurate  measurement 
of  the  sound  amplification  or  absorption  difficult. 

A  situation  more  amenable  to  experimental  study  results  if  the 
pumping  process  can  be  terminated  and  a  metastable  state  allowed  to 
develop  in  the  gas  before  the  sound  wave  is  introduced.  In  this  case  it 
is  obvious  that  the  decay  time  must  be  long  compared  to  the  sound  period. 
Three  typical  situations  where  such  a  condition  could  be  produced  by 
pumping  energy  into  a  particular  vibrational  mode  are: 


(1)  A  single  long  lived  relaxing  mode.  Most  stable  diatomic  gases 
fall  into  this  category.  Bauer  and  Bass  have  considered  the  case  of  a 
mixture  of  CO  and  H2  and  predict  large  amplification  when  the  period  of 
the  wave  is  close  to  the  relaxation  time  of  the  excited  state7.  As  men¬ 
tioned  above,  if  the  measurements  are  made  after  the  pumping  mechanism 
is  removed,  frequencies  in  this  range  are  ruled  out.  This  study  consi¬ 
ders  the  possibility  of  observing  the  amplification  when  the  sound  has  a 
period  that  is  short  compared  to  the  relaxation  time. 

(2)  Fast  v-v  and  slow  v-t  exchange.  For  most  pure  polyatomic  gases, 
and  in  some  gas  mixtures,  vibrational  energy  exchanges  rapidly  between 
modes  by  a  fast  v-v  exchange  and  then  more  slowly  with  translation  by  a 
slow  v-t  exchange.  It  was  early  recognized  that  such  a  "fast  series  pro¬ 
cess"  would  lead  to  single-relaxation-time  sound  absorption  and  dispersion 
curves  and  the  process  was  invoked  to  explain  the  single-relaxation  times 
observed  for  most  polyatomic  gases.8’9 

If  an  upper  vibrational  level  of  such  a  gas  is  pumped,  a  metastable 
state  with  different  vibrational  modes  at  different  temperatures  results. 
The  whole  manifold  of  vibrational  states  then  decays  at  a  slow  rate.  The 
resonant  v-v  exchange  within  a  given  mode  is  still  fast  compared  to  v-v 
exchange  between  modes.  An  example  of  this  stiuation  is  CH3F. 10,11  The 
vibrational  modes  equilibrate  on  a  time  scale  corresponding  to  100  gas 
kinetic  collisions  and  then  the  whole  system  decays  together  with  a  relax¬ 
ation  time  of  the  order  of  2  x  10~6  sec  atm  which  corresponds  to  about 
1500  gas  kinetic  collisions. 

Sound  propagating  through  a  gas  in  this  kind  of  metastable  state 
could  show  some  interesting  effects.  Even  if  the  sound  frequency  is  far 
above  the  v-t  relaxation  regions  and  is  not  communicating  translational 
energy  into  and  out  of  vibration  by  the  v-t  process,  there  still  may  be 
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communication  of  translational  energy  into  and  out  of  vibration  due  to 
the  non-resonant  exchange  between  vibrational  modes.  To  understand  this, 
consider  what  happens  in  the  compression-rarif action  cycle  of  the  sound 
wave.  When  the  translational  temperature  increases  upon  compression 
the  v-v  transition  rates  change.  The  number  of  molecules  traveling  the 
endothermic  direction  of  a  reaction  increases  more  than  the  number  going 
the  opposite  direction  with  the  result  that  one  vibrational  mode  increases 
in  temperature  more  than  the  other.  The  reverse  change  takes  place  in 
the  rarif action  part  of  the  cycle.  When  the  sound  frequency  is  close 
to  the  v-v  transfer  rate,  the  process  transfers  energy  from  the  compres¬ 
sion  to  the  rarif action  producing  sound  absorption.  When  there  is  an 
overpopulation  of  vibrational  states,  the  slowly  decaying  vibrational 
states  again  dump  energy  selectively  in  the  compression  part  of  the  sound 
cycle,  due  to  the  temperature  dependence  of  the  transition  rates.  This 
reduces  the  absorption  or  even  causes  amplification.  Measurement  of 
the  effect  could  give  valuable  information  about  relaxation  paths,  the 
transition  rates,  and  their  temperature  dependence. 

(3)  Slow  v-v  and  fast  v-t .  In  some  cases  the  v-v  rate  may  be  slower 
than  the  v-t  rate.  SO2  is  an  example  of  this  situation.12  The  v-t  rate 
is  about  10  times  the  v-v  rate.  In  this  case,  which  might  be  compared  to 
the  two  step  decay  of  a  radioactive  isotope,  energy  pumped  into  the  long- 
lived  mode  feeds  by  v-v  transfer  a  second  mode  that  is  communicating  with 
translation  by  v-t  transfer  and  maintains  its  vibrational  temperature 
above  the  equilibrium  value  for  periods  comparable  to  the  v-v  relaxation 
time  rather  than  the  v-t  time.  In  this  case  also  experimental  observa¬ 
tion  of  sound  amplification  may  be  possible  and  measurements  may  give 
valuable  knowledge  about  the  energy  transfer  process. 
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The  Reaction  Equation 

The  sound  absorption  (or  amplification)  and  velocity  can  be  obtained 
from  the  imaginary  and  real  parts  of  the  propagation  constant  for  plane 
waves  (ks).  For  this  purpose  we  write 


where  oj  is  the  angular  frequency,  M  is  molecular  weight,  R  is  the  gas  con¬ 
stant,  T  is  absolute  temperature,  and  Cv  and  Cp  are  the  specific  heat  at 
constant  volume  and  pressure.  In  the  presence  of  relaxation,  Cv  and  Cp 
are  time  dependent  or,  for  a  sinusoidal  wave,  complex  and  frequency  depen¬ 
dent.  The  sound  absorption  (as)  and  velocity  (vs)  are  given  by 

as  = Imaginary  part  of  ks 

vs  =  a)  x  [Real  part  of  ks]_1  2. 

The  problem  of  calculating  the  sound  absorption  and  velocity  is  thus 
reduced  to  finding  the  real  and  imaginary  parts  of  the  complex,  frequency- 
dependent  specific  heats,  Cv  and  Cp . 

For  this  purpose,  we  begin  with  the  reaction  equation  which  gives 
the  time  rate  of  change  of  energy  in  a  vibrational  mode.  We  will  assume 
that  the  normal  modes  of  vibration  can  be  treated  as  harmonic  oscillators 
and  that  resonant  exchange  of  energy  between  the  vibrational  levels  within 
a  single  mode  rapidly  establish  a  Boltzmann  distribution  within  the  mode 
characterized  by  a  "vibrational  temperature".13’14  In  the  metastable 
state  here  considered,  the  vibrational  temperatures  for  the  different  modes 
may  differ  from  each  other  and  from  the  translational  temperature  and  will 
be  changing  slowly  on  the  time  scale  of  the  sound-wave  period. 

The  time  rate  of  change  of  energy  in  the  a—  vibrational  mode  can  be 
1  5 


written  as 
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ia  -  ihujpjoglo  (k“^+1-M+J(a,b) 

-ka+i_>a,e+^e(a,b)  na+i(a)n6+i  (b) )  > 


3. 


This  equation  assumes  that  mode  "a"  is  exchanging  energy  with  a  second 
mode  Mb”,  and  with  translation.  Modes  a  and  b  can  be  in  the  same  mole¬ 
cule  or,  in  the  case  of  a  gas  mixture,  in  different  molecules.  The 


inclusion  of  additional  modes  of  interaction  is  straight  forward. 


l  5 


and  3  refer  to  the  vibrational  quantum  number  of  modes  a  and  b ,  h  is 

Planck's  constant,  Ua  is  the  fundamental  vibrational  frequency  of  mode 

/p  . 

a,  u  is  pressure,  is  the  total  number  of  molecules  carrying  mode  b, 
na  and  ng  are  the  number  of  molecules  with  a  and  3  quanta  in  modes  a  and  b, 
respectively,  ka  ^+^(a,b)  is  the  composition-averaged  rate  con¬ 

stant  in  the  dimensions  (time  x  pressure)  1  for  the  particular  process 
in  which  a  changes  to  a+i  and  3^3+j  in  the  mixture.  The  summations  over 
a  and  3  gives  a  net  rate  constant  for  transitions  where  mode  a  gains  i 
quanta  and  mode  b  gains  j  quanta.  In  general,  only  transitions  where  i 
and  j  have  different  signs  are  likely  to  occur,  i.e.,  mode  a  gains  i 
quanta  and  mode  b  loses  j  quanta  or  vice  versa. 

For  linear  perturbations  of  harmonic  oscillators, 

ka  a  ^(a,b)  =  {  [a(a-l)  •  •  •  (a-i+1)  ]  x  [  (3+j  )  (3+j-l) 

•••(B+l)]/i!j!}Xki4°’o4j(a,b)  4. 

In  terms  of  collision  efficiencies,  k3 >a+1 » p can  written 

k  =  EpXpMpP  (a,b,p),  5. 

where  P^  +ha»t>>p)  is  the  collision  efficiency  for  the  designated 

transitions,  p  denotes  the  collision  partner  (if  the  collision  does  not 
allow  the  transition  P  =  0) ,  and  Mp  is  the  collision  rate  of  a  molecule 
of  type  a  in  a  gas  of  type  p  at  unit  pressure.  MpPa^+*,^"^+ha»b,p)  is, 
therefore,  the  usual  rate  constant.  xp  is  the  mole  fraction  of  gas  of 


type  p. 
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We  now  assume  that  resonant  exchanges  of  energy  between  levels  of  a 
single  mode  rapidly  establish  a  Boltzmann  distribution  of  energy  within 
each  mode  so  that  we  can  write 

na  =  n0(a)e~ahVa/kTa  =  Na ( l-e-hVa/kTa) e"ahVa/kTa .  6. 

Thus,  each  mode  has  its  own  temperature  which  in  the  general  case  is 
different  from  the  translational  temperature.  Following  conventional 
notation  we  set 

=  0a,  and  l-e_9a/T  =  Xa(T)  .  7. 

k 


(Note  that  this  is  a  change  in  notation  from  ref.  15.) 

Using  equations  4  and  6,  the  reaction  equation  (Eq .  3)  can  be  summed 
over  a  and  g.  For  this  purpose  note  that 

£a(a+i)(a+i-l)-..(a+l)xa  =  ^  (-j-j)  =  (fi~) i+1 


After  summing  over  a  and  g  Eq.  3  becomes 

Ea  =  i?jiNahua  (?[k1+°5°^(a,b)e'10a/Ta  - 

[Xa(Ta)]1  [Xb(Tb)]j. 


o+i,  j-K> 


(a ,b) e 


-j9b/Tb 


]/ 


8. 


For  a  discussion  of  the  case  where  the  modes  are  degenerate,  see  ref.  15. 


A  Single  Long-Lived  Relaxing  Mode 

If  the  vibrational  temperatures  can  be  determined,  Eq.  8  allows  the 
determination  of  the  sound  absorption  and  velocity  in  the  excited  gas. 

The  three  cases  discussed  in  the  introduction  would  be  delineated  by  the 
choice  of  transition  rates.  To  illustrate,  we  specialized  to  the  simplest 
case,  that  of  a  gas  with  a  single  long-lived  relaxing  mode.  In  this  case 


j  0  and  i  1.  The  sum  over  i  and  j  has  only  one  term  and  Eq.  8  becomes, 
after  some  manipulation. 
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Ea  =  -  k10Xa(T) [Ea(Ta)  -  Ea (T)  ]  ,  9. 

where  Ea(Ta)  =  Nak0ae_0a/Ta/xa(Ta)  and  Ea(T)  =  Nak0ae~9a/T/Xa(T) .  This 
of  course,  is  the  expected  relaxation  equation  with  the  relaxation  time, 
t,  equal  to  [  k1 0 (Xa (T) ] -1 . 


Now  consider  the  case  where  a  sound  wave  is  propagated  through  a  slowly 
relaxing  gas.  Either  by  differentiating  Eq.  9  or  by  differentiating  Eq.  3 
and  the  summing  over  ct  and  g,  we  obtain 


E°  {  [i  ^  ^  .  AEa(Ta)  -  AEa(T) 


rEg 


} 


10. 


Lp  dT  k10  dT  T^Xa(T) 

Here  we  make  use  of  the  fact  that  the  transition  rate  is  a  function  of 
translational  temperature  only  and  we  neglect  higher  order  terms  in  the 
sound  variables.  Ea,  Ta  and  T  now  represent  unperturbed  values,  i.e., 
values  without  the  sound  wave.  Introducing  the  sinusoidal  time  dependence 
and  setting  AEa(Ta)  =  Ca(Ta)ATa,  AEa(T)  =  Ca(T) AT,  and  Eg  - - -Ea(Ta) ~Ea(T)  , 
this  equation  becomes 


(1+jorr) 


AT 


AT 


a  


=  [ 


Egx 


1  dp 


Ca(Ta)]  X  [P  dT  +  k17 


dk1 0 (a) 

dT 


]  +  -ga(T)  u 

TzXa(T) J  Ca (Ta)  11 • 


In  this  form  it  is  easy  to  see  the  effect  of  the  overpopulation  of  excited 


states  on  the  usual  absorption  due  to  relaxation.  In  the  "usual  case", 
that  is  when  the  gas  has  not  been  "pumped",  Ea  =  0,  and  Ta  =  T,  and  the 
right  hand  side  of  Eq.  5  is  unity.  Then  ATa/AT  =  (l+jorr)""1.  Writing  the 
frequency  independent  part  of  the  specific  heat  as  Cco,  the  frequency 
dependent  specific  heats  become 


Cv(w)  -  Coo+CaATa/AT  =  C00*fCa  ( 1+j ojt )  1 

and 


12. 


CpCm)  —  Cv(m)+R 

The  familiar  absorption  and  velocity  versus  frequency 
those  corresponding  to  a  single-relaxation-time  process. 


curves  are 


When  Ta  ^  T,  Eg  ^  0  and  the  frequency  dependent  specific  heats  are 


Cv(w)  =  4-  C^/(l+ja)T) 


0ae 


-0a/T 


with  =  E§x[-  +  d[lnk10(a)]/dT  -  +  Ca(T) 


and 


Cp  (co)  —  Cy(m)  +  R  — 


Egx 


Td+jcox)  ’ 
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1  dp 


Here  R  is  the  universal  gas  constant  and  ~  has  been  given  its  ideal  gas 
value  of  The  extra  term  in  Cp  results  from  the  fact  that,  at  constant 

pressure,  9p/9T  =  0. 

It  is  possible,  therefore,  to  plot  Cp(oo)  and  C v(w)  as  or>0,  and  the 
sound  absorption  and  dispersion  in  a  particular  gas  as  a  function  of  the 
decay  rate,  Eg .  Bauer  and  Bass  have  done  this  for  a  CO/H2  mixture,  treat¬ 
ing  the  CO  as  having  two  vibrational  states.  If  Eg  is  set  equal  to  the 
negative  of  their  pumping  rate,  the  two  treatments  give  similar  results. 

For  the  purposes  considered  here,  it  is  more  useful  to  set 


Ea  =  -r"l[Ea(Ta)  -  Ea(T)] 
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and  express  the  results  in  terms  of  the  vibrational  temperature  Ta.  Sound 
amplification  will  occur  when  Ca  becomes  negative.  (Note  that  Eg  is  nega¬ 
tive.)  We  now  seek  to  determine  values  of  Ta  needed  to  produce  this  sound 
amplification,  the  magnitude  of  the  amplification,  and  the  experimental 
conditions  under  which  it  can  be  observed. 


Calculation  and  Observation  of  the  Sound  Amplification 
For  most  diatomic  gases,  (the  halogens  would  be  an  exception),  the 
relaxation  time  is  so  long  in  the  pure  gas  that  collisions  with  impurity 
molecules  usually  control  the  relaxation  process.  For  such  collisions, 
^^dT10) iS  °f  the  order  of  a  few  one  hundreds16  and  amplification  will 


occur  when  Ta-T  is  a  few  tens  of  degrees. 
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However,  observing  this  amplification  in  a  metastable  state  in  the  absence 
of  pumping  is  difficult.  For  frequencies  where  the  sound  period  is  short 
compared  to  the  relaxation  time  (a  requirement  for  the  metastable  state  to 
persist  over  a  number  of  sound  vibrations)  ,  got  will  be  large  and  the  gain 
per  unit  length  of  sound  path  will  be  small.  Furthermore,  the  length  of 
the  sound  path  is  restricted  since  the  transit  time  for  the  sound  wave 
through  the  gas  must  be  small  compared  to  the  relaxation  time. 

In  the  case  of  polyatomic  gases  with  molecular  dipole  moments  and 
correspondingly  high  optical  absorption  coefficients,  it  may  be  possible 
to  extend  the  sound  path  length  by  pumping  the  gas  ahead  of  a  sound  pulse 
with  a  laser  beam  swept  through  the  gas.  However,  such  a  configuration 
does  not  appear  practical  for  the  diatomic  case  here  considered.  In  the 
first  place  to  absorb  the  laser  light  the  diatomic  gas  would  need  a  molec¬ 
ular  dipole  moment  and  in  addition  the  long  relaxation  times  would  require 
large  volumes  and  hence  large  power. 

This  means  that  high  vibrational  temperatures  are  needed  to  produce 
observable  amplification  of  the  sound  in  the  diatomic  case.  Table  1  gives 
values  of  the  gain  in  db/wavelength  for  some  typical  diatomic  gases  as  cal¬ 
culated  from  Eqs,  1,2,9,10  and  11.  In  these  calculations  got  is  set  equal 
100  and,  in  this  range,  the  gain  is  approximately  inversely  proportional 
to  cot .  However,  values  much  smaller  than  100  violate  the  metastable  state 
requirement  since  at  this  value  t  corresponds  to  approximately  16  sound 
periods.  The  value  of  —  lnkio  has  been  taken  to  be  0.01,  which,  as 
mentioned  above,  seems  to  be  representative. 
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Table  I.  Sound  amplification  in  db  per  wavelength  for  typical  diatomic 

gases  at  300°  K  with  elevated  vibrational  temperatures,  Ta.  The 
derivative  of  the  natural  logarithm  of  the  transition  rate, 
d ( ] nki o ) / dT^  has  been  set  equal  to  0.01.  The  amplification  is 
approximately  inversely  proportional  to  o)T . 


gas 

v  (cm1) 

Amplification  (db/A)  with  cot  =  100 

Ta  =  400°K 

600°K 

1000°K 

1500°K 

2000°K 

4000°K 

n2 

2331 

8.6  x  10-4 

0.016 

0.16 

0.54 

1.0 

3.4 

CO 

2143 

1.5  x  10-3 

0.024 

0.20 

0.61 

1.1 

3.6 

02 

1556 

7.5  x  10-3 

0.071 

0.36 

0.87 

1.5 

4.0 

C12 

557 

5.0  x  10~2 

0.26 

0.74 

1.4 

2.0 

4.6 
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Task  IV 

This  task  involves  general  studies  in  acoustic  cavitation.  We  have  made 
progress  along  four  lines: 

1.  In  a  study  of  the  nonlinear  oscillations  of  individual  gas  bubbles 
that  are  acoustically  levitated  in  a  liquid,  we  have  been  able  to 
measure  (a)  the  polytropic  exponent  of  the  gas  contained  within  the 
interior  of  the  bubble,  and  (b)  certain  aspects  of  the  pulsation 
amplitude  relating  to  excitation  of  the  second  harmonic  resonance  of 
the  bubble  oscillation.  Two  papers  describing  the  results  of  this 
particular  study  are  scheduled  for  publication  in  early  1983  in 

the  Journal  of  the  Acoustical  Society  of  America.  Some  of  this  work 
has  been  done  in  collaboration  with  Audrea  Prosperetti  of  the  Uni¬ 
versity  of  Milan,  Italy.  Abstracts  of  these  papefs  are  included 
As  Appendix  IV-A  and  IV-B. 

2.  Previously,  the  equations  that  describe  the  growth  of  air  bubbles  by 
rectified  diffusion  have  been  done  mostly  in  a  piecemeal  fashion.  We 
have  generalized  these  equations  by  extending  their  range  of  applicabi¬ 
lity.  A  paper  describing  this  work  has  appeared  in  the  November 
issue  of  the  Journal  of  the  Acoustical  Society  of  America.  This 
abstract  is  included  as  appendix  IV-C. 

3.  We  have  examined  the  effect  of  dilute  polymer  additives  on  the  acoustic 
cavitation  threshold  of  water.  We  discovered  that  the  addition  of  very 
small  amounts  of  these  additives  increase  the  cavitation  threshold  -  a 
desirable  effect.  Our  explanation  for  why  this  occurs  involves  the 

way  the  reduced  liquid  surface  tension  (the  additives  are  surface-active) 
affect  the  nucleation  conditions.  A  paper  describing  the  results  of  this 
study  has  been  submitted  for  publication  in  the  Journal  of  Fluids  Engineering. 


54 


The  abstract  of  this  paper  is  included  as  appendix  IV-D. 

4.  A  tremendous  amount  of  interest  has  recently  been  generated  by  certain 

suggestions  that  ultrasonic  scanners  used  as  diagnostic  devices  may  cause 
transient  cavitation.  A  recent  paper  by  Flynn  [J.  Acoust.  Soc.  Am.  72 , 
1926-1932  (1982)]  has  predicted  the  presence  of  acoustic  cavitation  by 
these  devices.  We  have  attempted  to  verify  this  prediction  experimentally. 
We  have  found  that  cavitation  can  be  produced  by  single  acoustic  pulses; 
however,  our  pulses  have  been  at  high  kilohertz  frequencies  rather  than 
the  megahertz  frequencies  associated  with  most  diagnostic  scanners.  We 
are  redesigning  our  apparatus  to  examine  these  higher  frequencies.  A 
report  of  our  active  work  in  this  area  was  recently  given  at  the  Bureau 
of  Radiological  Health;  an  abstract  of  that  presentation  is  included  as 
appendix  IV-E. 

Further,  a  general  review  article  on  acoustic  cavitation  was  presented  at 
the  IEEE  International  Symposium  on  Sonics  and  Ultrasonics.  This  review  paper  will 
be  published  in  the  proceedings  of  that  symposium.  An  abstract  of  that  presenta¬ 


tion  is  included  as  appendix  IV-F. 
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Appendix  IV-A 

The  Polytropic  Exponent  of  Gas  Contained  Within  Air  Bubbles  Pulsating  in 
a  Liquid.  Lawrence  A.  Crura,  Physical  Acoustics  Research  Laboratory, 
Department  of  Physics,  University  of  Mississippi,  Oxford,  MS  38677 


Gas  bubbles  that  are  driven  into  oscillation  within  a  liquid  pulsate 
with  amplitudes  that  are  determined  largely  by  the  thermodynamic  behavior 
of  the  gas  contained  within  the  bubble.  Theoretical  approaches  to  the 
solution  of  the  motion  of  the  radius  of  the  bubble  normally  involve  the 
poly tropic  exponent  of  the  gas.  Measurements  are  presented  of  this  expo¬ 
nent  for  three  different  gases  and  for  a  wide  range  of  values  of  the 
exponent.  Comparison  with  available  theories  indicates  good  agreement 
between  measurements  and  predictions  provided  the  gas  bubble  is  not 
driven  near  one  of  its  harmonic  resonances. 


To  be  published  in  the  Journal  of  the  Acoustical  Society  of  America, 
early,  1983. 
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Appendix  [V-B 

Nonlinear  Oscillations  of  Gas  Bubbles  in  Liquids:  An  Interpretation  of 
Some  Experimental  Results.  Lawrence  A.  Crum,  Physical  Acoustics  Research 
Laboratory,  Department  of  Physics,  University  of  Mississippi,  Oxford,  MS 
38677  and  Andrea  Prosperetti,  Istituto  di  Fisica,  Universita  di  Milano, 

I taly . 


Measurements  are  presented  in  this  paper  of  the  pulsation  amplitude 
of  an  individual  air  bubble  that  is  levitated  in  a  glycerine-water  mixture 
by  a  stationary  acoustic  wave  operating  at  a  frequency  of  22.2  kHz.  Obser¬ 
vations  of  the  bubble  pulsation  for  a  wide  range  of  bubble  sizes  demonstrate 

the  existence  of  the  n  —  2  harmonic  resonance  (to  z  to  /2)  .  The  available 

o 

theoretical  information  on  linear  and  nonlinear  bubble  oscillations  is 
adapted  to  apply  to  the  specific  experimental  conditions.  Comparisons 
made  with  the  theory  show  excellent  agreement  between  measurements  and 
predicti ons . 


To  be  published  in  the  Journal  of  the  Acoustical  Society  of  America,  early 


1983. 
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Appendix  IV-C 

Generalized  Equations  for  Rectified  Diffusion.  Lawrence  A.  Crum 
and  Gary  M.  Hansen,  Department  of  Physics  and  Astronomy,  University  of 
Mississippi,  Oxford,  MS  38677 

There  are  at  least  four  different  equations  present  in  the  litera¬ 
ture  for  the  threshold  for  growth  of  gas  bubbles  by  rectified  diffusion. 
These  various  equations  have  limited  regions  of  applicability  and  also 
give  significantly  different  predictions  in  common  regions  of  applica¬ 
bility.  In  this  paper,  a  set  of  general  equations  are  obtained  for  the 
growth  rate  and  the  threshold  for  rectified  diffusion  that  apply  over  a 
wider  range  of  interest.  It  is  shown  that  under  certain  restrictions 
the  previous  equations  agree  with  the  more  general  ones  given  here. 
Finally,  graphical  representations  of  the  equations  are  given  showing 

the  degree  of  approximation  and  the  range  of  applicability  of  the  spe- 
cialized  equations. 

Published  in  the  November,  1982,  issue  of  the  Journal  of  the  Acoustical 
Society  of  America. 
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Appendix  IV-D 


Effect  of  Dilute  Polymer  Additives  on 
the  Acoustic  Cavitation  Threshold  of  Water 


Lawrence  A.  Crum 
Professor  of  Physics 

and 

James  E.  Brosey 
Graduate  Assistant 

Department  of  Physics 
University  of  Mississippi 
Oxford,  MS  38677 


ABSTRACT 

Measurements  are  presented  of  the  variation  of  the  acoustic  cavitation 
threshold  of  water  with  concentration  of  the  polymer  additives  polyethylene 
oxide  and  guar  gum.  It  was  found  that  small  amounts  of  these  additives 
could  significantly  increase  the  cavitation  threshold.  A  theoretical  model, 
based  upon  nucleation  of  a  gas  bubble  from  a  Harvey-type  crevice  in  a  mote 
or  solid  particle,  is  developed  that  gives  good  agreement  with  the  measure¬ 
ments.  The  applicability  of  this  approach  to  an  explanation  of  cavitation 
index  reduction  in  flow-generated  or  confined  jet  cavitation,  when  polymer 
additives  are  introduced,  is  discussed. 


Submitted  for  publication  in  the  Journal  of  Fluids  Engineering. 
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Appendix  IV-E 

Acoustic  Cavitation  and  Some  Possible 
Biophysical  Effects 

Various  aspects  of  acoustic  cavitation  will  be  discussed  and  movies  will 
be  shown  of  the  different  types  of  cavitation  observed  in  a  liquid.  Recent 
observations  of  cavitation  effects  in  tissue  will  be  described  as  well  as  some 
possible  explanations  for  these  effects  in  terms  of  bubble  growth  by  rectified 
diffusion  and  transient  cavitation  inception.  Finally,  some  recent  work  on 
the  cavitation  produced  by  microsecond-length  pulses  will  be  discussed. 


Presented  at  the  Food  and  Drug  Administration  National  Center  for  Devices  and 
Radiological  Health,  Thursday,  December  9,  1982. 
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Appendix  IV- F 


ACOUSTIC  CAVITATION  IN  LIQUIDS  * 

LAWRENCE  A.  CRUM,  Physical  Acoustics  Research  Group,  Department  of  Physics, 
University  of  Mississippi,  Oxford,  MS  38677, 

A  broad  review  will  be  given  of  the  phenomenon  of  acoustic  cavitation,  generally 
defined  as  any  observable  activity  Involving  a  bubble  or  population  of  bubbles 
stimulated  into  motion  by  an  acoustic  field.  Discussions  will  be  given  of 
cavitation  inception  -  thresholds,  nucleation,  tensile  strength  of  liquids; 
bubble  dynamics  -  rectified  diffusion,  nonlinear  oscillations ,  cavity  collapse; 
and  cavitation  effects  -  sound  attenuation,  material  erosion  and  biological 
damage.  Movies  depicting  various  aspects  of  cavitation  phenomena  will  be  shown. 


*This  work  has  been  supported  in  part  by  ONR,  NSF  and  NATO. 


L.  A.  Crum 

Physical  Acoustics  Research  Group 
Department  of  Physics  &  Astronomy 
The  University  of  Mississippi 
Oxford,  MS  38655 
(601)  232-5840 


To  be  published  in  the  Proceedings  of  the  IEEE  Symposium  on  Sonics  and 


Ultrasonics . 
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Task  V 

This  task  involves  the  acoustic  measurement  of  energy  and  momentum 
accomodation  coefficients.  Under  this  task,  an  acoustic  technique  is  used 
to  determine  these  coefficients.  A  technical  report  was  distributed.  The 
abstract  of  this  report  is  included  as  Appendix  V-A. 
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Appendix  V-A 

Abstract  for  technical  report  submitted  October  20,  1982.  This  report 
will  be  published  as  a  paper  in  the  Journal  of  Chemical  Physics  in  early 
1983. 

ABSTRACT 

Acoustical  measurements  of  energy  and  tangential  momentum  accomodation 
coefficients  indicate  stable  and  reproducible  Pt  and  Ag  surfaces  are  obtained 
by  heating  in  a  vacuum.  EAC  and  TMAC  values  are  given  for  He  and  Ne  on  the 
heat-cleaned  Pt  and  Ag  surfaces  and  for  these  gases  on  Pt  and  Ag  surfaces  after 
exposure  to  various  adsorbent  gases.  AC  values  are  also  given  for  the  various 
adsorbent  gases  on  themselves  on  the  heat-cleaned  Pt  and  Ag  surfaces. 
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